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The two dimensional Heisenberg model for SrCu2(BOa)2 
has the exact dimer ground state which was proven by Shastry 
and Sutherland almost twenty years ago. The critical value 
of the quantum phase transition from the dimer state to the 
Neel ordered state is determined. Analysis of the experimen- 
tal data shows that SrCu2(BOs)2 has the dimer ground state 
but is close to the transition point, which leads to the unusual 
temperature dependence of the susceptibility. Almost local- 
ized nature of the triplet excitations explains the plateaus 
observed in the magnetization curve. 
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The pseudo spin-gap behaviors observed in high T c 
cuprates have stimulated intensive investigations on mag- 
netic systems with spin gaps. As products of this type 
of activities, several new spin gap systems have been 
found experimentally. Among them, some of the exam- 
ples which have two dimensional character include the 
coupled spin ladder systems, SrCu203 M, CaV2C>5 ||] 
and the plaquette RVB system, CaV^g || . 

Recently Kageyama et al found a new two dimensional 
spin gap system SrCu2 (1303)2 0- The crystal structure 
of SrCu2 (1303)2 is tetragonal and all Cu 2+ ions with a 
localized spin 5 = 1/2 are located at crystallographically 
equivalent sites. The two dimensional layers containing 
the Cu 2+ ions are separated by planes of Sr 2+ ions. 

The susceptibility drops sharply below the maximum 
at around T = 20 K. The peak of the measured sus- 
ceptibility is much suppressed compared with that of the 
theoretical value expected for a dimer model. The spin 
gap estimated by the nuclear magnetic relaxation rate is 
A = 30 K. Another novel feature appears in the mag- 
netization under high magnetic fields. The authors re- 
port that the two magnetization plateaus corresponding 
to 1/4 and 1/8 of the full Cu moment are observed for the 
first time for the two-dimensional quantum spin systems. 

The magnetic properties of SrCu2 (BOs)2 may be de- 
scribed by the two-dimensional Heisenberg model with 
the nearest-neighbor and next-nearest-neighbor cou- 
plings: 



h = j5>. Sj + J' ]T 



(1) 



The system is illustrated in Fig.^(a) and topologically 
equivalent to the model considered by Shastry and 
Sutherland || . We note that the nearest-neighbor bonds 
define a unique covering of all spins. On the other hand, 



the system with only next-nearest-neighbor couplings is 
equivalent to the square lattice Heisenberg model. We 
assume that the both coupling constants are antiferro- 
magnetic, J and J > 0, and then the present Heisenberg 
model is frustrated. 



• • 



• • 



1 S3 I 



I 



I 



• • 



• \ m - 



I 



I 



I 



• ♦ J m I • • 

(a) 



1 



J' 



4 



a 



(b) 



FIG. 1. 
SrCu 2 (B0 3 )2. 



(a) Lattice structure of the Cu 2+ spins of 
The nearest-neighbor bonds are expressed by 
solid lines and the next-nearest-neighbor bonds by broken 
lines. Square unit cells containing 4, 8, 16, 20 spins are shown 
by dashed lines, (b) Elementary unit for the interaction be- 
tween a pair of nearest-neighbor bonds. 

Shastry and Sutherland have shown that the singlet 
dimer state is an exact eigenstate of the Hamiltonian. In 
this paper we determine the critical value for the tran- 
sition from the gapful dimer ground state to the anti- 
ferromagnetically ordered gapless state. Analysis of ex- 
perimental data shows that SrCu2(B03)2 is a spin gap 
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system with the dimer ground state but is close to this 
quantum transition point. The unusual magnetic proper- 
ties of SrCu2 (B03)2 are explained consistently from this 
point of view. 

In the following analysis we use the dimer bases defined 
for each nearest-neighbor bond: 
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As Shastry and Sutherland 
product of the singlet states 



i*>=ni 



have shown, the direct 



(6) 



where a represents the nearest-neighbor bonds, is always 
an eigenstate of the Hamiltonian ([[]). 

The proof for the exact eigenstate is simple. Since the 
wave function is an eigenstate, actually the ground state, 
of the first term of the Hamiltonian, let us consider the 
effect of the second term. It is easy to show by elementary 
calculations that for any neighboring pair of the nearest- 
neighbor bonds the matrix element of the second term of 
Eq.(Q) vanishes: 



H ab\s)a\s) 



0. 



(7) 



To be explicit, Tt a b — J (si ■ S3 + S2 • S3) where the site 
indices are shown in Fig. [Tib) . Note that the vanishing of 
the matrix element is due to the odd parity of the singlet 
with respect to the reflection which exchanges the two 
spins Si <-> s 2 . 

The energy of the dimer state is given by -Edimer = 
— (3/8)7VJ where N is the total number of spins which 
is assumed to be even. It is clear that the dimer state is 
the ground state for small J <C J ■ 

Now we consider the other limit, J ^> J. In this 
limit the model is topologically equivalent to the two di- 
mensional square lattice Heisenberg model as mentioned 
before. According to the recent quantum Monte Carlo 
simulations the ground state energy per site is — 0.669J 
B. The first order correction due to the dimer coupling 
J is obtained from the spin-spin correlation of the next- 
nearest-neighbor pair of the square lattice model, which 
is calculated as 0.204 0. By using these results the 
ground state energy of the ordered phase is estimated 
as E AF = -7V(0.669J - 0.102 J). The transition point 
between the two phases may be obtained by equating the 
two energies, which leads to (J /J) c = 0.71. Since the 
ground state energy for the antiferromagnetic phase is es- 
timated by a variational calculation, the transition point 
thus obtained gives an upper bound. 

Figure ^| shows the energy obtained by the exact diago- 
nalization for finite systems, N = 8, N = 16 and N = 20. 



Periodic boundary conditions are used for the calcula- 
tions. By considering the finite size effects, we obtain an 
estimation, (J / J) c > 0.69. Therefore we conclude that 
the transition point is (J / J) c = 0.7 ± 0.01. 
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FIG. 2. Ground state energy per site for finite lattices: 
N = 8, 16, 20. The lowest ener are 
also shown by filled symbols. 

One can estimate the excitation gap in the dimer phase 
by the perturbation theory. After some calculations the 
spin gap up to the fourth order correction is given by 



2 V J ' 




(8) 



Up to this order the triplet excitations are completely lo- 
calized. This unusual behavior is understood by consid- 
ering the matrix elements for the triplet excitations. To 
be explicit we start from the state with a triplet with the 
spin quantum number s z = 1 at the bond a in Fig.|l|(b). 
The states connected with this state by the Hamiltonian 



n ab \h) a \ S ) b = J — |ti) |*o) 6 - — |t > a |*l>6 • (9) 

It is important to note that when a triplet moves to 
neighboring bonds it leaves another triplet behind. Next 
crucial observation is that the following matrix elements 
vanish also by the symmetry reason : parity with respect 
to the reflection, 



H ab\s) a \t m ) b = Q (m = 0,±l). 



(10) 



The above two facts, Eq.(g) and (|l0|), set a stringent 
constraint for motion of a triplet. Hopping of a triplet is 
possible through a closed path of dimer bonds and thus 
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the hopping processes start from the sixth order in the 
perturbation. On the other hand, in a similar system 
in one dimension, the triplet excitations are completely 
localized Q . In the present system, the triplet excitations 
are nearly localized with extremely small dispersion. 

The spin gap for finite systems is shown in Fig|| For 
the dimer region, J / J < 0.7, the finite size effects are 
small, which is a consequence of the almost localized wave 
functions of the triplet excitations. The perturbation re- 
sult given by Eq.(||) is very accurate for J / J < 0.5. 
On the other hand for the antiferromagnetic region the 
finite size effects are significant, indicating usual disper- 
sive magnon excitations. Since the dimer ground state 
is always an eigenstate it makes a level crossing with the 
Neel ordered state at the transition point. For such a 
case there are two possibilities concerning the nature of 
the transition ^ . From the results shown in Fig. 3 it is 
difficult to determine uniquely whether it is weakly first 
order or continuous. 




FIG. 3. Spin gap for fimte^SeS.Hf jyaftr&^tigl. The 

solid line is the perturbation result up to the fourth order. 
The dashed line is a fit obtained by adding fifth and sixth 
order terms to the perturbation result. 

Susceptibility at high temperatures is obtained by the 
high-temperature series expansion as 



4T 



J + 4J -J 2 + 8JJ +8 J 
4T + 16T2 



(11) 



From the expansion, the paramagnetic Weiss constant is 
given by 9 = (J + 4/)/4. 

Kageyama et al determined the excitation gap A = 
30 K ^from the NMR relaxation rate. It should be noted 
that from the low temperature increase of the suscepti- 



bility it is estimated as A = 19 K. In the following anal- 
ysis of the experimental data we use the former value 
which is identified as the spin gap by the authors of [0. 
The susceptibility data at high temperatures are fitted 
by a Curie- Weiss law with the Weiss constant 9 — 92.5 K 
and the effective <?-factor g — 2.14. The spin gap and 
the Weiss constant are sufficient to determine the cou- 
pling constants uniquely. By using the fit up to the sixth 
order shown in Fig.^, J = 100 K and J = 68 K are 
obtained. Temperature dependence of susceptibility is 
shown in Fig.^. In the figure comparison is made between 
the theoretical calculations for finite clusters N = 8 and 
N = 16 and the experiments. The theoretical values are 
calculated by the quantum transfer matrix method for 
the systems with periodic boundary conditions |Tc|| . Con- 
sidering the ambiguity for the estimation of the gap from 
the experiments on one side, and the smallness of clus- 
ters used for theoretical calculations on the other side, 
the agreement between the experiment and the theory is 
satisfactory. A better fit may be obtained when slightly 
smaller values are used for the Weiss constant and the 
spin gap. 
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FIG. 4. Temperature dependence of susceptibility for fi- 
nite clusters, N = 8 and lfj!9-i S=H i #§ , W%nd j = 68 K are 
used. The experimental data are shown by the thick solid 
line. Also shown are the theoretical results for smaller values 
of the ratio J / J = O(dimer), 0.2, 0.4 with the Weiss constant 
fixed 6 = 92.5 K. 

The characteristic feature of the susceptibility of this 
compound, the usual Curie- Weiss type behavior at higher 
temperatures above the peak and the steep drop be- 
low the peak, is well reproduced. The estimated cou- 
pling constants J / J = 0.68 is close to the critical value 
(J / J) c — 0.7. In Fig.|| also shown are the theoretical re- 
sults for smaller values of the ratio J / J = 0, 0.2, 0.4 with 
the Weiss constant fixed, 9 = 92.5 K. We can conclude 
that the closeness to the transition point is the origin of 
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the unusual temperature dependence of the susceptibil- 
ity. 

Next we discuss the magnetization curve. Figure |5| 
shows magnetization as a function of applied magnetic 
field. The results are obtained by the numerical exact di- 
agonalization for the system with N = 20 spins. One can 
clearly identify plateaus corresponding to 1/4 and 1/2 of 
the full moment. Experimentally, plateaus are observed 
for 1/4 and 1/8 of the full moment. Since the calculations 
were done for a small cluster, it is not possible to find a 
plateau at 1/8. The critical value for the plateau at 1/4 
estimated for the present set of parameters J = 100 K, 
J = 68 K is about 40 T which is comparable with the 
experimentally observed critical field H C 4 = 37T M. 




FIG. 5. 
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Formation of the plateaus are related with the almost 
localized wave functions of the triplet excitations. At spe- 
cial values of magnetization where the triplet excitations 
can take a regular lattice structure, energy cost would 
be locally minimum. It is natural to assume that the 
commensurability energy is more favorable when a unit 
cell is a simple square since the original structure has 
the tetragonal symmetry. The square unit cells are pos- 
sible for N = 4, 8, 16, 20, 32, • • • which corresponds to 
plateaus at 1/2, 1/4, 1/8, 1/10, 1/16, see Fig.|(a). 
The plateaus observed experimentally, 1/4 and 1/8, are 
just two of them. Our theory predicts that there will 
be a clear plateau for 1/2 of the full moment at around 
100 T and also for 1/10 at a smaller magnetic field than 
H C 2 = 28T. Since the plateau of the latter is expected 
to be small, a lower temperature and a clean sample will 
be required. The experiments on these plateaus are a 
crucial test for the present theory. Another prediction of 
the theory is the superstructure at each plateau which 



should be observable, for example by NMR. 

In conclusion we have identified that SrCu 2 (603)2 
has the exact dimer ground state found by Shastry and 
Sutherland ||. In the model the quantum phase tran- 
sition occurs at J / J ~ 0.7 from the dimer state to the 
antiferromagnetically ordered state. This transition is ei- 
ther weakly first order or continuous and SrCu2 (603)2 
is a unique spin gap system which is close to the quan- 
tum transition point. Unusual temperature dependence 
of the susceptibility is a consequence of the closeness to 
the transition point. Another interesting feature appears 
in the excitations. The low lying triplet excitations are 
almost localized and easy to form regular lattices under 
certain magnetic fields. The commensurability energy 
associated with the superstructures leads to the plateaus 
in magnetization curve at 1/2, 1/4, 1/8, 1/10, 1/16, ••• 
of the full moment. 
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